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Abstract. Numerical simulations of the evolution of strange-mode instabil- 
ities into the non-linear regime have been performed for a wide range of stellar 
parameters for Wolf-Rayet stars. It has been shown that the Wolf-Rayet models 
reach radial velocities which amount to up to 30% of their escape velocity. The 
acoustic luminosities suggest a connection to the observed mass loss. Most of 
the models show a jump in the mean effective temperature after reaching the 
non-linear regime. This jump is related to the run of the opacity. 

To determine the fate of unstable stars, the evolution of their instabil- 
i ties has to be fol lowed into the non-linear regime by numerical simulations 
( Grott et al. I [20051 ). Linear stability analyses have shown that strong instabili- 



ties with dynamical growth rates prevail in many types of massive stars . These 
are associated with so called strange modes ( Gautschv Glatzel Ill990l ). Mod- 



els of Wolf-Rayet stars exhibit th e strongest strange mode instabilities found 
so far ( Glatzel Sz Kiriakidis I fl993l ) . Thus they are promising candidates for a 



non-linear stability analysis. 

The Wolf-Rayet models investigated reach final radial velocity amplitudes 
in the non-linear regime of the order of 10 cm/ s. These values are comparable 
with the escape velocity. To investigate a possible connection between stellar 
instabilities and a pulsationally driven mass loss, the acoustic energy flux is 
considered. It represen ts the mechanical e nergy transferred by shock waves to 
the star's atmosphere ( Grott et al. 20051 ). We emphasize that the acoustic 



luminosity and thus any mass loss rate based on it strongly depends on the 
artificial viscosity which has to be used in the simulations. However, since 
viscosity has a dissipative effect the derived mass loss rates have to be regarded 
as lower limits (see Table 1). 



Table 1. Acoustic luminosities and mass loss rates based on them for mod- 
els with different initial effective temperatures. For the models, the mass- 
luminosity relation of the He-ZAMS was adopted. 
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After reaching the non-linear regime all unstable models exhibit a jump 
in the mean effective temperature. The final mean effective temperature is 
independent of mass and largely independent of the initial effective temperature. 
We find only two final mean effective temperatures around 30000-fT and 15000-fT 
respectively. Models with an initial effective temperature above « 32000.K" reach 
a mean effective temperature around 30000K, models with initial temperatures 
below this value end up with a mean effective temperature around f5000i"C. 
Considering the run of the opacity, these two temperatures correspond to the 
descending branches of the opacity from the respective maximum towards lower 
temperatures (see Fig.l). Should descending branches of the opacity in general 
give rise to a range of possible final mean temperatures of pulsating models, a 
third range above ~ 100000-ftT would be expected. In fact, some models show a 
jump to this temperature. The final mean position of the pulsating models in 
the HRD is of particular interest. Most of the models investigated move to the 
right and settle in a narrow temperature strip (Fig.l). It is remarkably close to 
the observed positions of Wolf-Rayet stars. 
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Figure 1. Left: Run of the opacity in the stellar envelope as a function of 
temperature. Different stellar models were used (different gray scales corre- 
spond to different models) to cover a wide range of temperatures. For further 
details see text. Right: Hertzsprung-Russell-Diagram. Indicated are the posi- 
tions of the initial models by crosses and asterisks and the final mean positions 
by squares and circles. The mean luminosities stay constant. Arrows indicate 
the direction of the evolution. For models with question marks, we were not 
able to determine a final mean temperature. Models without arrows don't 
change their mean temperature. The Hc-ZAMS is shown for comparison. 
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